The quantum effects in the undulator of infinite length 



V.G.Bagrov\ V.V.Belov^, M.M. Nikitin^, and A.Yu.Trifonov^ 

^High Current Electronics Institute, 

Siberian Division Russian Academy of Science 

4 Akademichesky Ave., 634055 Tomsk Russia 

'^Department of Applied Mathematics 

Moscow Institute of Electronic and Mathematics, 

B. Vusovsky 3/12, 109028 Moscow, Russia 

^Department of Mathematical Physics, 
Tomsk Polytechnical University, 30 Lenin Ave., 634034 Tomsk, Russia 

Abstract 

The first order quantum correction to the power of spontaneous radiation of electrons 
in an arbitrary two-component periodic magnetic field was obtained. The phenomenon 
of selfpolarization of the spin of electrons in a process of spontaneous radiation was also 
studied. By electron's motion in a spiral magnetic undulator, the quantitative charac- 
teristics of selfpolarization (the polarization degree and the relaxation time) are different 
from corresponding ones in synchrotron radiation. The limiting cases of near-axis and 
ultrarelativistic approximation were considered. 

Introduction 

The investigation of spontaneous radiation of electrons moving in periodic structures 
(undulators) is now an important and rather well developed branch of modern physics. A 
detailed bibliography about this problem can be found, for example, in reviews [1-2] and 
monographs [3-5]. 

In the majority of studies, the radiation was investigated by methods of classical elec- 
trodynamics [1-5]. Quantum-electrodynamical considerations were applied in few studies 
[1, 6-15], where the character of quantum corrections to the radiation of ultrarelativistic 
particles (7 ^> 1) was calculated and, in particular, it was shown that these corrections 
are small in actual undulators. 

Nevertheless, because of wide usage in theoretical and experimental research of mag- 
netic undulators, it is important to give a detailed analysis of the character of quantum 
corrections to the radiation over the whole range of electron energy allowed by the undu- 
lator regime. This range also includes the nonrelativistic energies, for which the effects 
caused by particles motion are significant. 

In the present paper we considered the first quantum correction to the characteristic 
of undulator radiation of charged spinor particles in a magnetic undulator. This consider- 
ation is limited to the case of a vector potential of an arbitrary magnetic field depending 
only on one spatial coordinate (along the particle drift). This particular model of an 
external magnetic field allows us to consider typical regimes of partical motion in mag- 
netic undulators: along the plane periodic trajectory for a plane undulator, and along 



the helical trajectory for a spiral undulator. On the other hand, this model allowes us 
effective calculate quantum-electrodynamical characteristics by using Bloch's stationary 
wave functions obtained in [16] by the one-dimensional WKB method [17]. 

Let us give some general conclusions. The results of our calculations [9-15] which 
are of certain theoretical interest for the problems of motion and radiation of particles in 
magnetic undulators of all energies when the quantum (spinor) properties are taken into 
account. 

1. If the characteristics of spontaneous radiation of electrons in a periodic magnetic 
field (calculated by methods of quantum theory in the first order over the radiation field 
when the motion is semiclassical) are represented in the form of formal (and generally 
asymptotic) Taylor series in the Planck constant h as h ^ 0, then each term of this series 
can be represented as a functional of the classical particle trajectory, the external field, and 
derivatives of this field on the classical trajectory. In particular, the n quantum correction 
contains the n-|-2-derivative of velocity, i.e., n+l-derivative of the vector potential. In the 
case of n = 1 [18], we calculated boson radiation in an arbitrary two-component periodical 
magnetic field, in which the helical motion is realized. 

It should be mentioned that the dependence of the characteristics of spontaneous ra- 
diation (calculated in the semi-classical approximation as h ^ with relativistic accuracy 
to the first order of h) on the parameters of classical trajectory was first shown in [19-20] 
for the case of electron energies in the relativistic range. 

2. The quantum corrections to the radiation characteristics and the infiuence of dif- 
ferent quantum effects (in the whole energy range allowed by the undulator regime of 
electron motion) depend^ on the specific form of the potential of the external field. For 
electrons moving along the same classical trajectory in different external fields, the quan- 
tum corrections (in comparison with the classical term) in general will differ. 

This conclusion is illustrated here by calculations of the characteristics of radiation of 
spinor particles moving in magnetic undulators. For example: 

1) The explicit dependence of quantum corrections to the frequency of the photon 
radiation on the specific form of potential of the external periodical magnetic field is 
shown (see Eq. (2.9)). 

2) The quantum expansion parameter of the total power of undulator radiation (in 
the on-axis approximation) qualitatively differs from the corresponding parameter in syn- 
chrotron radiation and depends explicitly on the character of magnetic field and its first 
derivative on the classical trajectory. It is interesting that, for a specific form of the field 
in which the particle moves along the spiral trajectory, in the ultrarelativistic limit the 
first quantum correction to the full radiation power coincides with the expression of the 
quantum term in the power of synchrotron radiation of the charge moving along the spiral 
in a constant and homogeneous magnetic field [21]. 

3) The quantitative characteristics of the effect of radiational self-polarization of elec- 
trons in a helical undulator (the degree of polarization and the relaxation time) differ from 
those in the theory of synchrotron radiation [22] and also from similar results in a field of 
a fiat electromagnetic wave [23] and in a axisymmetric focusing electric field [24], though 
in all considered (examined) cases a electron movement on a spiral is realized. This fact 
proves that the structure of the external field but not the classical electron trajectory plays 
a significant role in the phenomena of self-polarization. It should be mentioned that, in 
the ultrarelativistic limit (7 — > 0), the characteristics of the process of self-polarization 
calculated in Sec. 6 coincide (with an accuracy of ~ 7^) with those obtained earlier in [7] 
on the basis of an operator semi-classical method of calculation. 



In the ultrarelativistic case (7 — > 00), the dependence on the structure of an external field can be neglected. 
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1 Semi-classical electron wave function 



Let us consider an electron moving in a stationary magnetic field with the vector potential 

A{z + l)=A{z), ^ • ' 

where Ai{z) are arbitrary smooth periodic functions, I is the field period. Let us describe 
the electron motion on the basis of the relativistic Dirac equation: 

nD'^ = E^, Hd = c{a,V) + p^ruQC^ , (L2) 

V = {-ihV) - -A{z), 

c 

a \ / 1 



a VO -1 

where a = (ci, cr2, ^s) are Pauli matrices, E is the electron energy and e = — eo is the 
electron charge. Let us separate the solutions of equation (1.2) over a polarization states 
by condition: 

c St^ = Q\m^ St = {^',V), (L3) 
a 



a 



Here the parameter ^ describes the spin orientation with reference to the direction of 
motion (longitudinal polarization) either along the direction of motion ((^ = +1) or in the 
opposite direction (^ = — 1) [25]. 

Taking into account of quantum integrals of motion 

pi^=pi'$, P2^=P2^, p = —ihV, 

the solution of system (1.2) and (1.3) will be represented by: 

*£;,pj_,C = ^exp|^(pj_,r)|0£;^p.^_(^(z,;i), (1.4) 

where 

4'E,p^,({z,h) = {gE,p^,(:{z,h),XE,p^,(:{z,h)y, p±_ = (pi,p2,o), 

I k\ CcA , . ^ , . (1.5) 

XE,pU^z,n) = — KgE,p^,cy,z,n), r = {x,y,z), 

h + niQC^ 

N is a normalization constant. Spinor gE.p^,c{^^T^) satisfies the condition 

|('?,^±) + 0-3^3} 5£;,pj_,c(^'^) = C^9E,p^,dz,h), (1.6) 
^± = (7^1,7^2,0), 

V^=p^ + —Ai{z), i = l,2 
c 

or in the equivalent form gE,pj_,ci^^ ^) = {9+{z, ^, E,p±, (),g^{z, h, E,p±, ("))* 

{(P3 + CA)(Pl-iP2)~Hp3-CA)+Pl+ip2}<7+(^,^) =0, . . 

g.{z,h) = {Vi-iV2)-\C\-h)g+{z,h). ^ 



^In case this couses no misanderstanding the indeces z, ft, E,p^,C, may be omitted. 
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In the absence of turning points 

the WKB-solution of system (1.6) has the form [26]: 

g+{z,h) = f+{z,h)ex.p^^S{z)Y 



(1.8) 



(1.9) 



where f^{z,h) = f^\z,h) + hf^\z,fi) is a regular series over h —>■ 0. Here, the actual 
measureless expansion parameter is the ratio of the de'Brojlie electron wave length to the 
period of changing of field A ^ where [19] 



Tojo = 27r, c/?||r = /, 
Substituting (1.9) into (1.7) we obtain: 



dz 



nn E 

T = — 

Jo P{z 



z 

S{z) = J p{z)dz, 



(1.10) 



/f (z) = exp{- / dz 



/+\z)=if'^\z) I dz 



(0), 



(Pi - iV2) d {p{z) - CA) 
2p{z) dz {Vi - 1P2) 

{Vi+iV2) d ff\z) 



2p{z)ff{z)dz{ri-ir2) 



Here, the dot defines the differentiation with respect to the parameter z {C, = ^). By 
condition (1.8), the spinor g{z,h) = {g+{z, h), g-{z, h)Y , where 



, r(^iTj^)(ATCp)\V2 ri [ ( n 



CA {V2V1 - V1V2] 



{r±,r±) 



)dz}()l. 



11^ 



X < 



l-'l ( P ^iT^^2PTCA j {V2V1 - ViV2f + aV 



4 \p2 Vi^iV2 P 



2p%V±,V±) 



-dz 



H+ = l, M- = C, P = P{z) 
satisfies (1.6) with accuracy up to 0{1^). Prom here and from the condition 



9E,p^,c{^ + l,h) = exp|-g(£;,p_L,C)/| c/£,px,c(^'^) 
we find the quasi-momentum of an electron: 

q(E,p±,n = - / ip + C-^ — TT^ — fdz. (1.12) 

ijv ^2p {v±,r±) 8p^ {v±,r±) J 

In the absence of turning points p^{z) > 0, the normalization condition on the electron 
charge for the semiclassical stationary solutions (1.4) of the Dirac equation 

J ^+^bd^x = Sb',b, b={E,p±,C) 
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is equivalent to the condition 



AL^NN' 



1 + 



-L 



{E + moc^){E' + tuqc^) 
By analogy with the case without spin [12], we obtain 

L I 



d{q — q' 



-L 



Sqq' 



q-q'=0 



S 



E,E', 



d{E - E>) 

9E,p^,(9E,p^,-C = 0{h). 

Here, we used the fact that the quasi-momentum is defined with the accuracy up to ^j^n: 

L 



2E 



E + moc^ 



7 / dz g^^^^^^gE,p^,c 



d{q-q') 



-1 



\q-q'=0. 



d{q - q') 



d{E - E') 



q-q'=0 



1 [ 9E,p^,ci^)9E,p^,ci^)dz = j j'^\2X-Cn 



d{E - E' 

I r ^ / i — /-fi- vi'P2-'PiV2 \ 

V1 P2 - V1P2 \ 



As a result, for the normalizing factor we obtain 

2 ^ E + mpc^ 
4cA(2L)3 ■ 

2 The spectral-angular distribution of power 

The power radiated to the element of solid angle dO, = sin OdOdip at the electronic tran- 
sition from the state to the state with the radiation of photon can be found by 
usual electrodynamical methods and is equal to [21]: 



dWK, 



dn 



00 



(2.1) 



{e-K , Bhfi' ) , Oct = (co- , Bhfi' ) ; 



Bb,b' = J ^'^,a^bexp{-i>c{n,r)}dx, 

= (cos if cos 9, sin (p cos 9, — sin 9), 
e„ = (sin (f, — cos 9?, 0) , 
n= (61,62,63), 
ei = cos (^sin^, 62 = sin(^sin0, 63 = cos 0, 



(2.2) 



where a are the Dirac matrices and hxn is the momentum of the radiated photon. We 
assumed, that the average electron drift along the axes Ox and Oy is absent: 



/■^i(^), f 



p{z) 



dz = 0. 



(2.3) 
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The spectral-angular distribution of the total radiation power can be obtained from 
(2.1) by summation over all final states E',p'j_: 



dn 



E 



dWb,t 



By calculating the matrix elements (2.2) the wave functions (1.4) are represented in a 
form of Bloch functions and their periodical parts are expanded in Fourier series. After 
integration over coordinates, this leads to the following conservation laws: 



q — q' = eshx + 27rhn/l. 
For the spectral-angular distribution of radiation power we obtain: 



(2.4) 



r,2 OO 



n=l 



dn 



27r 



*(>f)=0 

I 



{|(e^, B{n, C, C'))l' + l(e., ^(n, C, C'))l'}^-5) 



,1 



S(n,C,C') = {2LfNN'- I dze-'-'^'4>%,,^^,^,{z)a4>E,p^,^{z), 



$(x) = h~^{q - q') - xes - 27rn/L 



(2.6) 



The radiation frequency u{n, h) can be defined from the system of conservation laws 
(2.4) and the conservation law of energy E — E' = chx. For this, we find the increment 
the of quasi-momentum Ag = q — q', which, according to (1.11), is equal to 



Aq = hwpoil) + ^(C - C')V3(0 + ^^Vpi(0 + \h^ujC,'p2{l), 



where 



X 



dz 

p{zy 



(2.7) 



(2.8) 



Pi{x) 

P2{x) 



cx 



E 



-{e±,P±)] -p'ei 



dz 
p^{z)' 



V2C1 - Vie2 V2V1 - V1V2 



E{r^,V±) ie±,V±) 



AV 



dz 

p{zy 



Pz{x) = - 



A f V2V1 - V1V2 



dz. 



(2.9) 



Taking into account the condition (2.3), we will get po(0 — Substituting (2.7) 

into (2.4), we find the radiation frequency (up to 0{h^)) in the form: 



u{n,h) = |i - f -^('^cipi(0 + C'p2(0)| ; 

a^cl = ^{c^n-^(C-C'>3(0}, 



(2.10) 
(2.11) 
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By analogy with (2.7), we obtain 



1 + 



^(^/>l(0 + |/92(0) 



(2.12) 



<I>(3<)=0 

It should be noted, that the expression for the frequency of classical radiation (2.9) 
thus obtained differs from the corresponding expression in [12] by a summand, which is 
proportional to (C — C')P3(0 ^.nd corresponds to transitions with spin flip. The proba- 
bility of these transitions is proportional to h. Therefore, at ^ = 0, the frequency (1.9) 
difi^ers from the frequency of classical radiation in the field (1.1) to harmonics whose the 
probability of radiation is equal to zero. This means that they just coincide. 



3 The matrix element of the transition currents 



Substituting (1.5) into (2.6), we obtain 



B{n, c, CO = ^(1 + CC')^^^(^, + ^(1 - CO^T^n, C), 

I 



B^Hn,0 = '-^-^J I dzexp{-ixe3z}g+,^^^_^{z)agE,p^,c{^ 



(3.1) 
(3.2) 

(3.3) 



Let us consider the matrix elements without spin flip 



gl{z,E',p^,Q g±{z,E,p^,C) \ _ 
gUz,E'J^,0 g^{z,E,p^,0 



f^^\z,E',i/^,Of'±'{z,E,p^,C) I 1 + 



f^^^*{z,E',f^,Of^^'{z,E,p^,0 I l^iCh 



f(0) 



hV2Vi-ViV2 
2 {rf + rl)p 



(P^CA) 



> X 



p2 ± 



xexp{ - / dz\p{z,E,pi_) - p{z,E' 



> ; 



/f*(z,E',pl,C)/f (z,^,p-l,C) 



p 2p y p- 



ie±,P±)- 



E 



^. eiPj-eiVi ^ E/c-{e±,V±) E_ 
^ Vl + rl p cA 



X ex.Y>{iC,hHzp2{z)/2\; 

f^^\z,E'J^,C)f^'\z,E,p^,0 = - 



X ^ ( ±C\p 



E (g±,P±) 
cA Vf + Vi 



Vi ± iV2 - ^/ixjei ± ie2 + (Pi ± iV2) 
-(ex,^±)))} 



E 
c 



X eyiTi>{iC,h>czp2{z)/2} 
and by analogy with the scalar case, we have 



z 

j \p{z, E,p±) - p{z, E'J^)]dz = iozpoiz) + ^hiuhpiiz). (3.4) 



7 



So, within an accuracy up to the quantum correction of the first order over h ^ 0, we 
obtain: 



' J P 





p. 



(3.5) 



-ei - CA 



iP2 / E {e±,P±) 



I f dz 



I J p 



V2 + — SV2 ^ 62- (3.6) 



1 f i„ 

ij 




>cp^ 

{ E/c- {e±,r±) 



p 



(3.7) 



where 



+1 63 + <A— 2-— 2— > 

HP vi + vi ], 



ipi{z, n, C, h) = ujzpo{z) - ujz^ + ^J^zpi{z) + C^L0zp2{z) 



(3.8) 



and pj(-2;) are defined in (2.7). The matrix elements (3.4)-(3.6) have the structure 



sTT(n,C)=s"(n,C)-^4 / 7 exp (po(^) - ^) } . (3.9) 



The last term does not influence the radiation power. So we will consider B^^nX) = 
B {n,Q. For the expressions with the spin flip, we perform analogous calculations: 

/f *(z,^',p-l,-C)/f (2,i?,pl,C) = ±{V±,V±)'^^expm2iz,C)-xe3z)}, 

P^' {z,E ,p^,-C)f^'{z,E,p±,C) = ±C (-p2^-p2)i/2 exp{z(V'2(^,C) - >^ze3)}, 

where 

ip2{z,C) =uj(^zpo{z) - ^z^ + C,zp3{z) = uFq{z) + C,zpz{z). 



(3.10) 



We obtain: 



I 



hx vriQC 1 f iV2^ — CPiP 



;1 f iP2^ 
iJ P{VI 



2 AMy pipi + viY/-" 

fiK uiqc 1 f —iXPi — Cp'P2 



2 I J pivi + viyi'^ 







2 A2 I 



I 



P 



eyi^{iip2{z,C,)}dz, (3.11) 
exp{#2(^,C)}(^^, (3.12) 
exp{z'02(2, C)}c?'2- (3.13) 



The obtained expressions (2.5), (2.8), (2.10), (3.1)-(3.11), in principle solve the prob- 
lem of taking into account the quantum corrections caused by quantum effects of motion 
itself hjjQ I {E (3\\^) , and by output of radiated photon fvuj/E. 
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4 The full radiation power in the on-axis approx- 
imation 



To integrate the spectral-angular distribution over all angles and to sum over the spec- 
trum in order to obtain the total radiation power with an accuracy to the first quantum 
correction is possible only in the near-axis approximation. It can be characterized by the 
classical parameter 



u = max 

te[o,T] 



(4.1) 



which is the maximum deflection angle of the electron velocity from the undulator axis. 
Here c/3j = (?VjlE is the velocity along the j is the axis, T is the period of motion over 
trajectories 



dz 



T=^r 

c2 J p{z) 





which is related to the period of field I hy I = c(3\^T, where c/3|| = ^ is an average 

drift along the axis Oz, and V = iVi-,V2,p) is a kinetic momentum. 

Summing the spectral-angular distribution of the radiation power over a spin and 
integrating over the time in (2.8), (3.1) and (3.4)-(3.6), we use the following relations: 



dz ^' dt 



1,3. 



Then 



E 

Vi = -A, 

c 



J F{z)dz = J F{z{t))cMt. 



--(e±,P±) = -(l-(/3±,e-l)), 



(4.2) 



zpq{z) - e-iz = ct- (fci(i),n) = Fo{t), zpi{z) = -^Pi{t); 











$(x)=0 





E 



nujQ 
-00 ' 



Wo 



2tt 



where 



Piit) 



{[l-{P±,e^)?-Plel}^, 



and for the matrix elements of transition currents 

T 

B{n) ^ 



fc = (0,0,1), ex = (ei,e2,0), P± = {Pi, p2,0), 
ei = sin 9 cos tp, 62 = sin 6 sin (p, 63 = cos 6. 



(4.4) 
(4.5) 
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The final result for the radiation power into the element of solid angle = sin 9 dOdip 
is given by 



(4.6) 



dfl 27rc w^i^n 

71=1 CI 

aa{n) = {B{n), e^), a^{n) = {B{n),eg), i/jq = 1 - /Jyea, (4.7) 



= (sin(^, — sin ip, 0), 
eV = (cos sin ^, sin(/7C0S ^, — sin^), 



where c/3|| is the average drift velocity along the axis Oz. The unit vector n = (ei, 62, ea) 
points in the direction of photon emission. 

The expressions (2.1), (4.2)-(4.4) obtained above solve the problem of including the 
first quantum correction into the radiation of a spinlcss particle. It should be noted, that 
the same result can be derived by using the solution of the Klein-Gordon equation in the 
field (1.1) [12]. 

Using the properties of conjugate Fourier series [27] in expression (4.5) it is possible 
to sum over n and to obtain the spectral-angular distribution of the radiation power 
considering the first quantum correction in the form: 



W 



^fdn 

27rc / 2T 



—^dt + ^ J dh J dt2A{ti)D{t2) cot ; 



Oin"), (4.8) 



S = 7r[c(ti - t2) - {(n,fo(ti)) - (n,ro(i2))}], 



where the double integral is understood in the sense of the principal value, and the func- 
tions D{t),A{t) have the form: 

A^{t) = {A{t),e^f + {A{t),ee)\ A{t) = [/3(n,/3') + V/3']/V''; 
A{ti)D{t2) = {A{ti),e^){D{t2),e,) + {A{ti),e^){D{t2),e^y, 
d B{t) 



dt V ' 



V' = F^ = l-(n,/3); 



B 



{t) = { (1 - (e"^, /3l))/33l(t)V' + [2/33^V + 2/33/3(n, P) - 0m 



63 

/3| 



+ 



m)9 d'z tx 1 



}: 



[l-(/3±,e^)]^ 

0i 



■dt 



WclV'O ' 



^ = p"^^ + p^{fi, P") + 3(n, + 0{n, 0')]. 

In the classical term of (4.6), one can integrate explicitly over angles, however, in the 
quantum term, this integration is impossible. 

Complete calculation of radiation characteristics can be done in the near-axis approx- 
imation /X — ^ (4.1). The terms of first order in |Lt — ^ in (4.2)-(4.5) imply 



/33 ~ /? « , /3i ~ A ~ , 
t 

Fo{t) = i^ot- J (ex, 0)dt + 0^(/i'), 



Pi(i) = /?n 



(l-Plei)t-2 ie^,P^)dt 
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t 

Fi{t) = -/?Vo j {e±J±)dt + O^ii?). 





Therefore, in the near-axis approximation we have 



9$ 



T 

B = — — / dt e^p{iuont}P 



0(.<)=O P\\ I ^0 J 

6 = na.ot-^(l + ^eo)+0, 



(4.9) 
(4.10) 
(4.11) 



V'o 



(i + =|sa)+o>^), 



where = huj/{Ef3j). 

For the spectral-angular distribution of the radiation power, we obtain 

n=l 

where 



5.(^,<^) = A(^) 
Saie,v) = faie) 



T 

^ J (A sin 99-/32 cos (/p)e^^«"*dt 


T 

^ y (/3i cos ip-p2 sin (^)e^'^°"*dt 



-cos^)2 



After integration over angles, we find the spectral distribution of radiated power i 
the following form: 



in 



,2, ,2 °° 



c 



'■K + So-) ^ 



n=l 



tJ 



5 _1 1 

'-'Tr — o 



+ 0^{^^^) + 0{h^)= (4.13) 

oo 2 ^ j. 

^n2(5, + 5,)^ - / Hk{z)e"--'dz 

n=l k=l Q 



+ 0^(/x2) + 0(n2), 



5.= 



3(1-/32)2 
1 

(1-/3J)2 



1 5 + 19/3j? 



(4.14) 



1 - ^0" 



1 + 3/32 



where -ff(^;) = (— ^2('2;)5 .4.1(2;), 0) is the magnetic field and a;„ = 2iTn/l. 

For summation over the spectrum we use the Parseval equality and the properties of 
conjugate Fourier series [27] were used: 

_ 9 _ 



2 00 

EE 

1=1 n=l 



n 







T 
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2 oo 

EE 

1=1 n=l 



n 



iujont 



Pidt 



1 1 



1 1 

JdtJ dr(/3i(t),/3l(r))cot(^7r) = ^<5i. 







The last integral is understood in the sense of the principle value. So, the total 
radiation power with account of the first quantum correction is 



^ n 5i 5 + 19/3g 
£;/3j<5o20(l -/Jj)' 

3 h 6i 3 + 9/3g 

4 6o 4(1 - pfA ' 



(4.15) 
(4.16) 



3c(l-/?2)2T 



3m2c3(l-/32) 



If a component of the magnetic field is equal to zero, then we have the radiation power 
in the plane undulator regime [9, 10]. 

Remark. From a ratio (4.11) it follows that in on-axis approximation the maximum in 
distribution on harmonics it is necessary that number n for which factor H(n) of expansion 
in a Fouier series for a magnetic field force is maximum. In first, on visible, this conclusion 
verifying the qualitative Motz results [28] (in the respect that not without fail maximum 
on the module is first (n = 1) Fouier factor) was made in work [29]. In particular, from 
(4.11) it follows that of magnetic fields, for which the influence of a harmonic H(ncr) is 
essential, where Ud (1 — P'^^)/£,o, the semiclassical expansion of (1.9) type is not true, 
as the quantum amendment becomes comparable with classical summand. 



5 The radiation power in the hehcal undulator 

In general, it is impossible to integrate the spectral angular distribution over all angles and 
to sum it over the whole spectrum. The total radiation power including the first quantum 
correction can be obtained in the ultrarelativistic case only for special magnetic field, 
in which the electron moves along the helical trajectory. Let us consider the periodical 
magnetic field 

H = { — Hq cos az, —HQsmaz,0} (5-1) 
with the period I = 27r/a. The relativistic equation of motion is 

x" = z'luq sin az, y" = —z'lvo cos az, (5-2) 
z" = uJo{x' sin az + y' cos az), 
luq = cqcHq/E, eo = |e|. 

In the general case (with the arbitrary initial conditions), the solution of system (5.2) 
is expressed by the elliptical integral. However, in the case of interest, when the condition 
of absence of the average drift along the axis Ox and Oy is fulfilled and the projection 
of initial momentum on this axis is zero, the dependence of coordinates on time is of the 
form 

f{t) = {Rsmuot,—RcosuJot,cP^i^t}. (5.3) 
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The system (5.3) describes the motion of a particle along the helix with frequency of 
spiral motion uq = 2tt/T = c(5\\a and the radius of spiral R = (5±/{(3\\a), [5± = eoH/(Ea). 

In the field (5.1), the spectral-angular distribution of radiation power is defined by the 
general formulae (4.4) and (4.5) and is of the form (see Appendix Al): 



W 



\a^{n)\ 



2 _ 



n=l n 



COS 6 — 



sin6'd6' 



cos 0)3 



[|a^(n)|2 + |a^(n)|2], 



(5.4) 



sm( 



;^n^W(i3^ + S)- (5-5) 



Z 

2 

7-2/ 



1 - 



K(n)|^ = /?x^7„^(z) 



z 

hcoon 

^2 



1 + 2^11 COS z 



/3|| cos 



+ -1- 



-|(i-y^n(^)j;(^ 



where 



n/9x sin 9 



d 



Jn{z) = -rMz), 



1 — /3||cos0' 

and Jni^) are the Bessel functions [30]. The expression (5.4) as ^ — coincides with the 
well-known formula for radiation power for the motion along the helix in a homogeneous 
magnetic field [21]. 

The integration over angles and the summation over all frequencies in the quantum 
term of expression (5.4) can be done only in the ultrarelativistic limit j3 1, j3'^ = f5^^+j3j_. 
We change to the system of reference, which moves along the initial axis Oz with velocity 

[21], and we use the asymptotic representation of the Bessel function Jn{z) and its 
derivatives as the argument z — > — [21]. The result for the total radiation power with 
account of the first quantum corrections is 



55^3 h f E 



2elc 



16 rriQcR \mQC^ 
E 



(5.6) 



3 \moc^ 

and the /^ri fa snoe n- and a-components of polarization which are equal to: 

E 



_ 1 5^/3 h 



8 16 mocR \moc^ 



_ 7 50\/3 h 



8 16 rriQcR Vmoc^ 



E 



(1-/3 ), 



(l-/3y)- 



6 The probability of transition with the spin flip 

The probability of transition w{C,, (') with the spin flip can be found by analogy with the 
expression (2.5): 



2n 



^^'^'^ = ^/^^ / ^ ^ x{|aTi(,,c)|^ + |4i(n,C)|^}, (6.1) 

n=[uo] 

aV-{n,0 = {e^,Bn{n,0),qquadaV-{n,0 = (e^, 5Ti(n, 0), 
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where uq is determined by tlie condition x(fo) = 0. In matrix elements B^^{n,Q it is 
convenient to change the integration variables dz = {(?p/E)dt. The result is 

V'2(t) =a;Fo(t) + CP3(t), i^o(0 = ^M-(rci(*)>^)]> (6-3) 



Now let us consider the model of a helical undulator (5.1), which is characterized by 
a helical electron trajectory: 

relit) = — {P± sina;oi, -P± coswot, u^oPut}. (6.4) 



In this case 



tp2(t) = moot - zsm{ip - uot), p3{t) = -^ujQt; 

P\\ 

cqH cqHI eoHl 

PI 



Ea 2TTE 2TTmQC?'^^ 



CP 

z = zo [n - — ] , zo 



I3±_ cos t 

and it is possible to integrate over the time [13] 

Mn, = - /?')^/'e--'^ {CPiiM^) + P-^Jn{z)] , (6.5) 

«a(n, C) = "21^(1 - /?')'/'e-*"^ { (/? j,(z) + CP\\ ^Mz)) cos 9 + +CP±Jn{z) sir(e.)) 

Integration of the expressions over angles and their summation over the spectrum is 

possible in the case of the near-axis approximation, which is characterized by the maximal 
deflection angle of electron velocity from the undulator axis n = 0\\ (4.1). In the case of 
helical undulator this means that the parameter zq = (/3_l sm.6)/ipQ is very small. Keeping 
in (6.4), (6.5) only terms, which are liniar over zq <C 1, we obtain 

/?~/33~/3||, = ^o{n - Q) / il^o; 

a,(n,C) = - /3^)V2e— |i/?(l - C)<5„,_i- 

-\0\\z5n,O + ^/3(1 + C)^n,l + \P{1 - Oz5n,2^ ; 

-(Z/3COS0- 2C^xsin0),5„,o+ [/3(1 +C)cos0 + C^/?isin0]5„,i + ^(3{l + Q cos 6 8^,2] ■ 
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Here, we used the following asymptotic representation of Bessel functions [30] : 
After summation over the spectrum we have 

t „2, ,3 

-(C,-C) = ^^|4/3i{A + M' (6-7) 







fa= r \,e ' sin^dg, (6.8) 



J ^0 



Finally, we have: 



^(C, -0 = on^^m„ n H^v^ 1 1 - , ^ ' (6-9) 



30^2 (1-/32)3 1^ ^5 _ 2/32 + 3/34 

C) = ^{1 - cr(/3 

15£;2/34(i-/32)3c 



^(C,-C) = ^{i-Cr(/3)}, (6.10) 



The function 



" hw^Plelib - 2/32 + 3/34) 



5 - 2/32 + 3/34 
is a monotonic function within the limits: 

1 = r(o) > r(/3) > r(i) = o. (6.11) 

It follows from (6.6) that the probability of the electron radiation with spin flip is a function 
of the initial orientation of the electron spin. This fact causes the effect of radiational 
selfpolarization of spin in a bunch of electrons [21, 22]. This means that, under certain 
condition 

«1, /?x = ^, /3^ = l-(^)^ (6.12) 



/3(1 - /32)V2 ' El' ^ E 

independently of the initial orientation of the spin of particles in a bunch, the primary 
spin orientation along the electrons motion will be determined. It follows from (6.7) that 
the polarization is absent, as /3 — ^ 1, but it is complete as /3 — ^ 0. However, the condition 
(6.9) does not allow us to consider these limiting cases. 

Let us consider in (6.6) the relativistic particles with /3 — > 1. Then, we find the 
relaxation time r and asymptotic degree of polarization p: 



These expressions within an accuracy up to 0(7~^),7 » 1 coincide with expressions 
previously obtained in [7] . 
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Appendix. The first quantum correction of the 
boson radiation in a hehcal undulator 

After the summation over the spin, the spectral-angular distribution of the radiation 
intensity of the electron moving in field (5.1) with the first quantum correction has the 
form 

1 n 





27r n(3± sin 9 

T 1 — Po cos t' 

where 

|a^(n)|2 = \Bi{n) cos e - Bsin) sine \aa{n)f = |S2(n)p; 
-B(n) = / (i?7 exp |i ( nr? — zi sin ^ ^o^t ^ ^.^ \ 1 ^ 



(Al) 



1 - TT^l S z(l - /?± sine cos ; 

[ 2EP^ 1 - Po cos 6* J 

We will integrate (A. 2) over rj using the following relations: 





27r 

i_ y gi(nr,-«sin,7) sin^ 7] di] = -Jn{i)] 


27r 


27r 


i_ y e*(nr,-«sin,) gin r? COS 7? dr? = i| - ^Jn(e)' 



27r 

^ j g«(n,-«sin,) si^2rysinryd7? = 2i + "^^MO - \jn{0 



2-K 

^ y g^(n,-€sin,) sin 2,7 Sin 77 dr/ = -2n |i - + ^Mo] 
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and 



(3 = {/3_L cos r], -/3± sin rj, Po} ■ 
As a result for the matrix elements we obtain 



n 



Bi{n) = I3± { —Jn{zi) ( 1 - ^ ^ ) + "1^1 o a 

\z\ \ 2EPq 1 - Po cos ) 1 — Po cos 



1 /n^ 1 • \ 

-P±sme i-^Jnizi) - -Jnizi)] + 

— 2ra^ ri^ + 1 • 1 ■■ 



i?2(n) = -iP^ 1 



n 



2S/?2 1 - /3o COS ""'^'^ ' m 1-/30 cos<? 
^/3iSine jn(2i) - ^Mzi)^ + 

+2^2 ( -Jnizi) - -^Jn{zi) + ^Jn{zx) 



z Z^ Z^ 



Bs{n) = /3o 1 



n 



2EP^ 1 - Po cos 6* 



?fcJO n 

"E;a^ l-/3ocos6» 



ft ' " f Til ' Tl \1 

P± sine Jn{zi) + 2Z2 ( -Jn(-Zl) - —Jn{zi))\ | • 



-2z • ■ \z 

Substituting (A.4) to (A.l) we obtain (5.4) 



(AA) 
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